Superconductivity and Ferromagnetism in Oxide Interface Structures: Possibility of 

Finite Momentum Pairing 
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We introduce a model to explain the observed ferromagnetism and superconductivity in LAO/STO 
oxide interface structures. Due to the polar catastrophe mechanism, 1/2 charge per unit cell is trans- 
ferred to the interface layer. We argue that this charge localizes and orders ferromagnetically via 
exchange with the conduction electrons. Ordinarily this ferromagnetism would destroy supercon- 
ductivity, but due to strong spin-orbit coupling near the interface, the magnetism and superconduc- 
tivity can coexist by forming an FFLO-type condensate of Cooper pairs at finite momentum, which 
is surprisingly robust in the presence of strong disorder. 



Introduction - It is known that a conducting electronic 
state can form at the interface between two insulating 
oxides [T]. A particularly well studied example is the Ti02 
interface between SrTiOa and LaAlOa. The carrier den- 
sity can be controlled by a backgate on the SrTiOa side 
and superconductivity (SC) has been discovered over a 
range of densities with maximum T^ of about 0.3 K[5]. 
Recently, signs of ferromagnetism (FM) have also been 
reported[SHZl- In particular, Li et aL|S| showed that SC 
and FM coexist in the same sample and that the FM 
moment is large, sa 0.3 — 0.4/xb per interface unit cell. 
Assuming that the FM and SC arise from the interface, 
these observations raise the question of whether the SC 
has to be unconventional in order to coexist with FM. Be- 
fore addressing this question we need to understand first 
the nature of the electronic state at the interface and up 
to now no clear picture has emerged [51 |H |TT] . Are most of 
the electrons localized or extended? Does the FM come 
from local moments or the mobile electrons and what is 
its origin? In this paper we propose a model for the inter- 
face electrons which is consistent with existing transport 
data. Based on this model we explain the existence of 
FM and the coexistence of SC and FM. For the latter, 
the key idea is that a large Rashba-typc spin-orbit cou- 
pling exists at the interface [T2]. Such a Rashba term 
is particularly favorable for the formation of a conden- 
sate at finite momentum, called a Fulde-Ferrell-Larkin- 
Ovchinikov (FFLO) state which coexists with FM p^ITi] . 
This general idea was pointed out earlier by Barzykin 
and Gorkov[T5]. However, they considered only the clean 
case and their solution is quickly destroyed by disorder. 
Surprisingly, with increasing disorder the FFLO state is 
revived P^. We suggest that the SC observed at the in- 
terface is described by this disordered stabilized helical 
FFLO state. This state is sometimes referred to as a 
"helical" FFLO stateTJ] since, pairing occurs at a sin- 
gle momentum q so that A(r) = Ae'^i"" unlike the usual 
FFLO state where pairing occurs at both ±q so that 
A(r) = Acosq • r. 

The Model - As shown in Fig. 1, LaAlOs consists 
of layers with alternating charged, while SrTiOa has 




O O O O O O (AlOi)' 

(LaO/ 

O O O O O O (AIO2)" 

(LaOf 

000000 (A102)" 

(LaOf 

Interface O -O*- O -&^ O -0>- O (TiOj)'" 

SrO 
z Ti Layer 1 O O O O • • • TiO, 



L^-xTi 



Layer 2 Q O O O O O 



SrO 
O TiOo 



FIG. 1. Schematic depiction of the SrTiOa/LaAlOa oxide 
interface structure. Filled and empty circles depict Ti and 
Al ions. The half electron charge per unit cell is transferred 
to the interface Ti02 layer localizes and orders magnetically 
(shown as arrows on the interface layer) via exchange polar- 
ization of conduction electrons on the subsequent Ti layers 
(shown as wavy cloud in Ti Layers 1 and 2). Inset: Disper- 
sion of electron bands arising from 3d orbitals on Ti layers 
near the interface. 



charged-neutral layers. As a result of the charge dis- 
continuity at the interface, an electric potential propor- 
tional to the number of LaAlOa layers is built up. This 
phenomenon is termed the polar catastrophe. Since the 
Ti ions allow for mixed valence charge compensation, to 
avoid the polar catastrophe, half an electron per unit cell 
is transferred from the surface AIO2 layer to the Ti02 
across the interface. The electrons are expected to oc- 
cupy the dxy orbital on the Ti atoms. Due to the rela- 
tively narrow bandwidth, an on-site repulsion U and a 
nearest-neighbor Coulomb repulsion V will cause these 
electrons to be localized on every other interface site. 
This picture of local moment formation at the interface 
has been proposed before [^. Super-exchange via the 
oxygen is expected to provide a weak antiferromagnetic 
exchange. 

The application of a back gate or the presence of de- 
fects forces more electrons to the interface. We assume 
that the effect of U and V makes it too costly to place 
these electrons at the interface layer. Instead, the ad- 
ditional electrons mainly occupy the Ti layer next to it 



(called layer 1), and their wavefunctions spill over to lay- 
ers 2 and 3 as well. We model these electrons as occu- 
pying two dimensional conducting bands. These are the 
electrons seen in transport measurements, with a typi- 
cal areal density of 3 x lO^'^cm"^, or about 10% of the 
density compared with the localized electrons. These 
electrons initially go into the dxy band [TU]. Hall ef- 
fect measurements show nonlinearity in magnetic field 
H with increasing gate voltage, which has been inter- 
preted as the appearance of a second carrier with lower 
mobility [ini HZl [TH]. We assume that these are the d^z 
and dyz bands, illustrated in the inset of Fig. 1. These 
bands are highly anisotropic with a heavy mass in one di- 
rection which may be responsible for the lower mobility. 
The dxz and dyz bands are higher in energy because their 
lobes point towards the negative charge at the interface 
and because their bandwidths are narrower. 

In addition, transport measurements show that the 
elastic scattering rate 1/r drops rapidly with increas- 
ing carrier density |18j . Furthermore, from the analysis 
of magnetoresistance, it was found that a Rashba term 
i^R = az-{(Txk) grows rapidly with the gate voltage [H], 
reaching a spin splitting Ago — 2akF of 10 meV near the 
peak of the superconducting dome, a value comparable 
to the Fermi energy of « 40 meV. Since a back gate volt- 
age tends to pull carriers from the interface, the decrease 
of 1/r is reasonable but the increase of Ago is counter- 
intuitive. We believe this trend is a consequence of in- 
creasing admixture of the dxz and dyz bands with increas- 
ing carrier density. The Rashba energy is determined by 
the polarization of the electron wave-function due to the 
asymmetric environment at the interface and the contri- 
butions come mainly from near the atomic core, where 
the electron is subject to a large electric field —dV/dz: 
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where the wavefunction tpk.. (r) — J^t '^e4>k (*") and £ de- 
note various angular momenta which are admixed due 
to the asymmetric environment of the interface [19j . Let 
us restrict ourselves to admixtures between d, p and s 
states. Since \/V ex z, non-vanishing contributions in- 
volving the d bands in Eq.(l) come only from the cross 
term between d^z and p^ and between dyz and Py. Fur- 
thermore, the dxy band can have a nonzero Ago only via 
the admixture of either s and p or, more importantly of 
dxz with Py and dyz with px ■ We expect that the latter 
hybridization between d and p orbitals gives rise to A^o 
that increases with chemical potential. 

The final ingredient of our model is the exchange cou- 
pling between the local moments and the conduction elec- 
trons. We write the standard phenomenological form 
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where s{r) = ^tpl^rTajSippir) is the electron spin density 
operator in the dxy band, and Si is the local spin operator 
on site i. We introduce a similar coupling J^ for the dxz 
and dyz bands. It is useful to introduce Jq = Jk/it-q 
where uq is the inverse of the interface unit cell area, and 
similarly Jq — J^/np. The Schrieffer- Wolff expression is 

Jo = 2P{jj^ ^) where i is the hybridization between 

the local moment and the conduction band orbital and 
Ed < is the orbital energy of the local moment relative 
to the chemical potential. We note that the dxz and dyz 
orbitals in layer 1 are orthogonal to the localized dxy 
orbital in the interface layer, so that the hybridization f 
vanishes except for the admixture of other orbitals in the 
dxz and dyz bands. We therefore expect Jq <^ Jq- 

The Origin of Ferromagnetism - We note that the 
present problem is the opposite limit to the familiar prob- 
lem of dilute Kondo impurities, where the Kondo screen- 
ing of the local moments competes with RKKY interac- 
tions between them. Here the density of local moments 
rii = i^uq is much greater than the conduction electron 
density rij,, i.e., the separation between local moments 
^/y/ni is much smaller than kp^. We can still view the 
local moments as interacting via RKKY interactions, but 
this interaction will be ferromagnetic and with a rela- 
tively long range of (2A:i?)~^. The FM ordering tempera- 
ture Tp can be worked out ^20j and apart from a numer- 
ical constant, the result was shown to be equivalent to a 
mean field treatment of Hj which we shall adopt below. 
In this picture, which was introduced by Zener |21j and 
referred to as the Zener kinetic exchange mechanism, the 
local moments order by polarizing the conduction elec- 
trons. This mechanism has been applied successfully to 
explain the FM of Mn substitution in GaAs and we bor- 
row the results here [52] . We introduce the average local- 
ized spin order per site S = ^ Si('^i) ^^^ the average 
electron density s = :^ J dr{s{r)). To quadratic order 
the total free energy density takes the form 
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The last term is the mean field decoupling of Eq.(2). In 

the first term Xo = 5 ^° T+e where /lo = gfJ-B and S — 
1/2, g = 2 in what follows, 9 > is the Weiss term due 
to the weak AF super-exchange exchange which we shall 
ignore below. In the second term Xc — |:/^o^(0) where 
^{O) = vn* j-nh^ is the density of states including spin 
of a free electron gas. (The presence of a Rashba term 
does not change the spin susceptibility of a free electron 
gas [13 US].) Minimizing Eq.(3) with respect to S leads 
to a purely quadratic term in |«Sp, and the sign change of 
its coefficient determines the FM transition temperature 
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We find that J^(0)/no 



12 " no JT-o 
= 0.64(m* /TOe) eV~^ For 



m%y/m = 0.7, Jo — 1.3 eV will give the observed 
Tp « 300 K. By comparison, for Mn/GaAs, Jq is w 1 eV. 
Here i is smaller, but \ed\ is also smaller because the same 
orbital is involved in the local moment and the conduc- 
tion electron, so the estimated Jq appears reasonable. 
Thus we conclude that the Zener kinetic exchange mech- 
anism can account for a robust FM state. 

Next we estimate the polarization of the conduction 
electron. In the mean field theory, the effect of S on the 
conduction electrons is described by an effective Zeeman 
field Hj ss J druoHuF • •s(r) where Hmf = Jo^S. Li 
et aZ.[S] reported an ordered moment of — O.SfiB per 
interface unit cell, i.e., 0.6/xb per local moment in our 
picture, which implies |«S| = 0.3. Using Jq = 1-3 eV, we 
estimate a Zeeman spin splitting \^qHmf\ ~ 200 meV, 
which is comparable to or exceeds the Fermi energy. Thus 
the dxy band is largely spin polarized. Since Jg <^ Jo, we 
expect a smaller (but still significant) polarization of the 
dxz and dyz bands. 

Nature of the Superconducting State - We assume 
that the superconductivity originates from a conventional 
electron-phonon coupling mechanism, which is modeled 
by an attractive short range interaction g with a cut-off 
given by the Debye frequency ujd- Since the supercon- 
ducting transition temperature is Tc ~ 0.3 K, the pairing 
gap ~ 0.04meV is the lowest energy scale in the problem. 
In the dxy band we estimate a Zeeman splitting of 0.2 eV, 
which exceeds the Pauli limit by more than 3 orders of 
magnitude, and precludes the possibility of pairing in the 
dxy band. The dxz and dyz bands will also be partially 
polarized due to the exchange interaction. However, the 
exchange splitting, /igS, in these bands is expected to be 
much smaller (although likely still /io-B 3> A). Moreover, 
as we argue above, we expect that the Rashba spin-orbit 
coupling, Ago, is even larger in the dxz and dyz bands 
than that observed in the dxy bands. It is natural to 
look to strong spin-orbit coupling to preserve pairing in 
the dxz and dyz bands despite large Zeeman splitting. 

The enhancement of B^ due to Rashba spin-orbit cou- 
pling was first demonstrated by [15j the case of weak 
Rashba coupling (Aso ^ £f) and no disorder, and later 
by [12] for the case of weak Rashba coupling and mod- 
erate disorder. They showed that an FFLO state is fa- 
vored, where the pairing occurs with a finite center of 
mass momentum |24) . Here we extend their analysis to 
treat arbitrarily strong values of A^o and disorder. We 
begin by neglecting disorder and find the susceptibility 
to form Cooper pairs at finite pair momentum q = qy. 
The dispersion for the ± Rashba branches is: e^, .^ = 

(k+q/2)~ 



Rashba bands and 



-l/k 



tan-i(^). The key is that 
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H ± ay/k^ + {ky + q/2 + /ioS/a)2. [25] In the 
physically relevant limit: vpq, B <C A^^, and we can ex- 
pand in q and B: 



'k+q/2 ~ ^k 



e^{B^O) + ivFq/2±^ioB)sm4>k+0{vFq^,B^) 



(5) 
where vp — \J o? + 2^/m is the Fermi-velocity for the 



-k+q/2 ^-k+q/2 



choosing q ~ ^" [26; ensures e.^ 

0{^2-) foi' 0,11 angles 4>k- This should be contrasted with 
the usual FFLO case without spin-orbit coupling, where 
the linear terms cannot be cancelled for all angles for 
any choice of q. However, we cannot prevent an energy 
mismatch in both bands simultaneously. By choosing 
q = 2fj,oBy/vF we optimize for the e_ branch, which has 
a larger density of states, v^, and find: 



;UoB, « Ao(A,o/Ao)^+°/"^ 



(6) 



where Ao = cj£) exp 



is the superconducting 



gap in the absence of the Zeeman field fJ-oB. In the limit 
ck ^ Wi? we recover the results of Barzykin and Gor'kov: 
fioBc ~ Ao\/Aso/Ao. For the oxide interface system we 
expect stronger spin-orbit coupling, A^o sa sp, and find 
an even larger enhancement: iiqBc ~ Ao(Aso/Ao)^/^. 

The above calculation is only valid in very clean sys- 
tems for which Aot ^1. In practice, we expect to be 
in the dirty limit, Aot <C 1. To incorporate impurity 
scattering, we consider spin-less, short-ranged impurities 
and compute the disorder averaged Cooper-channel sus- 
ceptibility in the e^^r 3> 1 limit, by summing the ladder 
diagrams for impurity scattering (called the Cooperon). 
As shown in Fig.[2[ there are 3 regimes. First, in the weak 
disorder regime (r^^ < Ao) we find that the critical field 

drops rapidly to the Pauli limit Be ~ Ao(Aot) i-^/^'-f . 
This can be understood as follows: as in the clean case, 
the pair momentum minimizes the effect of the magnetic 
field in the e_ branch, while keeping the pair breaking in 
the e_|_ branch. However, impurities can scatter Cooper 
pairs from the — band to the + band where they rapidly 
decohere. Thus, disorder enhances the dephasing effects 
of the Zeeman field, which becomes fully pair breaking as 
r^^ — > Ao. On the other hand, for very strong disorder, 
r^^ ^ Aso, the Rashba bands e^ lose their identifies 
due to the rapid impurity scattering. Here, spin and mo- 
mentum become decoupled and the problem reduces to 
that of conventional parabolic bands with effective spin- 
orbit scattering rate f^""^ = A^^t <^ t^^. This is the 
D'yakonov-Perel limit where the spin diffuses in small 
steps between rapid impurity scattering(27j. In this limit, 
it was demonstrated in [55] that f^oBc « Ao/v^Aqt^, 
and SC occurs at g = 0. 

The interesting limit is the intermediate regime Ag <C 
T^^ <C Aso [16]. Here the disorder is weak enough that 
the Rashba bands maintain their identity, but a pair in 
the e^ band can readily be scattered to a pair in the 
e^ band. Unlike the weak disorder case, the pairing is 
strongly admixed and it is not possible to sacrifice the 
pair breaking of one band in favor of the other. On the 
other hand, due to the spin-orbit coupling both spin (the 
Zeeman energy) and momentum show diffusive behavior. 
As a consequence, the pair breaking effect of the magnetic 




FIG. 2. Critical Zeeman splitting, Be, measured with respect 
to the bare SC gap Aq as a function of disorder strength 
T~^ . The FFLO state identified in 15 is rapidly destroyed as 
Aqt — > 1. Remarkably, for stronger disorder the FFLO state 
re-emerges and Be is enhanced beyond the Pauli limit (shown 
as a dashed line). 



field is weaker. The pairing interaction is dominated by 
the Cooperon, C, which ordinarily develops a diffusion 
pole C = ■ "L 2 with D = ^^ , but in the present case q 

is replaced by g± ^° for the ± bands respectively. Due 
to the strong mixing of the pairing channels, the effec- 
tive pair breaking strength is given by the combination: 
Sa=± i^ {"^pq + A2/X0-B) . Since zy+ 7^ i^- this combina- 



tion is minimized by finite momentum q 



4aB 



and 



we predict that an FFLO state exists in the intermedi- 
ate regime, but with a different q from the clean case. 
The corresponding Bf, is: fJ-oB^ ~ ^— - Physically, on 
this limit tb "^ t and, unlike the weak disorder, the pair 
dephasing time t'^/t grows with increasing disorder. 

Discussion - To summarize, we propose a model 
to explain the coexistence of SC and FM observed in 
STO/ALO oxide interface structures. In this model, 
a half-charge per unit cell is transferred to the inter- 
face layer and forms a lattice of local moments due to 
Coulomb repulsion. These local moments then order fer- 
romagnetically via exchange with lower-density bands of 
mobile electrons residing in Ti layers near the interface. 
The large FM exchange would ordinarily kill SC in these 
mobile bands. However, the presence of a large spin-orbit 
coupling enables the formation of an FFLO state which 
can coexist with strong magnetism. In this FFLO state. 
Cooper pairs form with finite pair momentum perpen- 
dicular to the direction of magnetic ordering. Unlike the 
usual FFLO state without spin-orbit coupling [121 HIj 
spin-orbit coupling parametrically enhances B^ beyond 
the Pauli limit and enables FFLO pairing to survive to 
much stronger disorder (up to A^qT > 1). Experimen- 
tally, T~^ decreases with back-gate voltage Vq while si- 
multaneously Ago increases. In our picture, the observed 
Tc vs Vq dome is related to the behavior of B^ shown in 
Fig. [2] sweeping from r~^ > A^o to r^^ < Ago- 

We wish to emphasize that our model suggests that 
the FM and SC occur in different bands, and are uni- 
formly distributed at the interface. Another possible ex- 
planation of the coexistence is in the spatial separation 



of the two phases. Indeed, a recent experiment [B] shows 
that the FM is arrange in domains and might not occupy 
the entire interface layer. However, the large total mag- 
netic moment of the system [S] , implies that the domains 
should occupy most of the area at the interface. The SC 
order parameter, in contrast, does not go to zero at any 
point in the plane. 

Before concluding, we briefly discuss some experimen- 
tal signatures of our model. First, if the exchange cou- 
pling for the dxz and dy^ bands is such that the conduc- 
tion electrons polarize in the opposite direction of the 
interface moments, then the external in-plane field Hh 
would align the local moments along H\\ but would reduce 
the total Zeeman field seen by thed^^^ and dyz bands to 
^o(-B — i?||). This leads to the unusual prediction that Tc 
should exhibit a maximum at finite Hh . Second, one can 
look for q ^ pairing by creating a Josephson tunnel- 
ing junction with a conventional SC film. By applying a 
magnetic field parallel to the junction. Cooper pairs tun- 
nel at finite momentum Ak, and the Josephson current 
would peak when Ak = q. 
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Supplementary Materials 

In the main text we have explained that the conduction electrons at the interface can either occupy the dxy band or 
the dxz and dyz bands. According to our model, the electrons in the d^y band are polarized by the localized electrons 
and participate in establishing the the Ferromagnetic order. The exchange interaction between the electrons in the 
dxz and dyz bands and both the localized electrons and those occupying the dxy band is much weaker. Then, these 
electrons can undergo a phase transition into the superconducting state. Still, the exchange field plays an important 
role as it create an effective Zeeman field parallel to the interface, which is expected to exceed the Pauli limit. As we 
argued before, the strong spin orbit interaction induced by the asymmetric environment at the interface, can stabilize 
the superconducting state at magnetic fields much higher than the Pauli limit. This is done by pairing electrons with 
finite momentum (the FFLO state). 

In this supplement we present a sketch of the derivation of the gap equation from which we determine the critical 
exchange field for the FFLO state. In addition to the effects of spin-orbit coupling and the exchange field, we wish 
to examine the effect of disorder on the state. The full details of the derivation and its extension to various special 
cases will be presented in a longer companion paper [1]. 

We start from the following Hamiltonian describing the conduction electrons in the dxz and dyz bands: 

H = Ho + i/inip + H,j + Hsc, (7) 

where Hq is the Hamiltonian of the electrons in the presence of Rashba spin-orbit interaction: 

— fc2 



Ho = }_^ — 4(k)ca(k) + az ■ {(Tap X k) cj,(k)c/3(k) 

k.a 



(8) 



The operator cj,(k) {ca(k)) is the creation (annihilation) of an electron with spin a along the z — direction and 
momentum k. Here we assume that the kinetic energy of the electrons can be approximated by a parabolic band. In 
reality, the dxz and dyz bands will have different mass in the x and y directions, however, this does not qualitatively 
change the subsequent analysis (as we demonstrate below). We consider a random impurity potential: 



H,^^^J2^ir)ci{r)cair). 



(9) 



which has only short-range correlations: V{r)V{r') — V^S{r — r'). The effective Zeeman field B = Hmf is chosen to 
be parallel to the x-axis. Thus, 



Hj = ^ /-toB • (TQ/3c]j(k)c^(k') 

k,Q,;3 



(10) 



Finally, the pairing interaction is: 

Hsc^-g J2 4(q/2 + k)cljq/2-k)c^(q/2 + k')c_Mq/2-k'). (11) 

k,k',q,a,/3 

For the next stages of the derivation it wiU be helpful to transform into the chiral basis that diagonalizes Hq -\- Hj\ 
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FIG. 3. Diagrammatic representation of Cooper pair propagator (top line) and Cooperon ladder (bottom line). Straight lines 
with arrows are electron Green's functions, wiggly lines represent attractive phonon mediated interactions, dashed lines with 
X represent disorder scattering, and Greek letters label the ± Rashba bands. 



where 'i/'i (V'±) ^-I'S electron creation (annihilation) operators in the chiral basis. Correspondingly, the eigenvalues of 
electrons in the two chiral bands (also referred to as the Rashba branches) are: 



ExO^) 
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+ XaJkl + [ky + ^oB/a), 



(13) 



with A = ±1. In the chiral basis, both the impurity scattering matrix element as well as the superconducting 
interaction acquire momentum dependence: 



yA,A'(k,k')- 
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i{ky + fioB/a)'^ kx + i{ky + ^^B/ay 
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H^c^-Y. .9Ai,A.,A3a4(k,k',q)^l^(q/2 + k)7/^l^(q/2-k)7/^A3(q/2 + k')^A,(q/2-k'), (Mb) 

k,k',q 

where repeated indices are implicitly summed and. 



3Ai,A2,A3,A4(k, k',q) = 



kx + qx/2-i{ky + qy/2 + ^qB/uY k^ - Qxfi - iiky ~ qy/2 - fioB/a) 



|k + q/2 + MoSy/a| 



+ A2 



\k^q/2-fXoBy/a\ 



. kx+qx/2 + i{ky + qy/2 + ^oB/ay k^ - qx/'^ + i{ky - qy/2 - ^iqB /af 
-^3 r, — ; 77;-, r,^ I I \- M 



|k + q/2 + ^oSi//a| 



|k - q/2 - ^l^By/a\ 



(14c) 



We are now equipped to study the vertex equation for the interaction in the Cooper channel. This equation takes 
the form: 

rAi,A2,A3,A4(k, k',q) = gAj_A2,A3,A4(k, k',q) + 2_^ 5Ai,A2,j)ii)2(k, ki, q)/C^i^,,2,»)3:>)4(ki, q)r^3,^4,A3,A4(ki, k', q). (15) 

ki£i 

Here /C denotes the two particle Green's function in the particle-particle channel. In the presence of impurities, K, 
must include scattering by the disorder potential. Therefore, in the ladder approximation, the two particle Green's 
function is a product of two Green's function connected through the Cooperon, C',^^,j2,,73,r;4(qj£n): 



/C, 



VI ,V2 .Va -.m 



(k, £„, q) = Gr,, (q/2 + k, e„)G^2 (q/2 - k, -£„)C»?i,»?2,»?3,»)4(q: ^«), 



(16) 



where 771^2 are not implicitly summed. The single electron Green's function can be written in terms of the eigenenergies 



of the two branches given in Eq. 13 

GA(k,e„) 



i£n ~ Ex{k)+ n+ — sign(e„) 



(17) 



The quasiparticle life time is a function of the density of states of both bands, 1/t = V'^'k{v+ + i/_) with v\ = 
m(l — Xa/vp)/2Tr. The vertex equation is illustrated diagrammatically in Fig. ^ The critical exchange field is 



determined by the largest field for which there is still an instability towards pairing, i.e. for which the effective 
interaction in the Cooper channel still diverges. For |q|,/zoi3/a ^ kp, we can neglect the dependence of the bare 



vertex (Eq. 14c I on q and B. Then, the gap equation takes the following simple form: 



Ai,A2,A3,A 



,(k,k',q) =gAi 



■^2,^31-^4 



(k,k',q) 



1-.9 



/C, 



ii-i-,vi,V2,n2 



(ki,q) 



ki.ei 



(18) 



Note that although the vertex only includes contributions where the two incoming (outgoing) electrons are from the 
same Rashba branch, the Cooperon contains sections in which the two electrons occupy different branches. We will 
not show here the general solution for the vertex equation, as it is rather cumbersome. Instead, we will present the 
expression for the vertex in four different regimes and extract the corresponding values of the critical exchange field 
at T = 0. We wish to point out that since the exchange field enters the energy by shifting the momentum ky inside 
the spin-orbit term (see Eq. 13), the FFLO state have finite momentum along the y direction, q_ = qy. 



PAIRING INTERACTION VERTEX WITHOUT DISORDER 

We first analyze the Cooper vertex in the absence of disorder (r^^ = 0). In this limit, the Cooperon is equal to 



unity, and the vertex described in Eq. 18 diverges when 



1-9^2 [G+(q/2 + ki, £„)G+(q/2 - ki, -e„) + G_(q/2 + ki, £„)G_(q/2 - ki, -£„)] = (19) 

ki,e„ 



One can see that the onset of pairing occurs simultaneously in both bands. To simplify the integration over the 
momentum and frequency, we expand the single particle energies in powers of q/kp and B/akp: G\{(\/2 + k, £„) w 
+ Aa|k| + /i + {vpq — A/igi?)]^^. To linear order, the interaction vertex becomes: 



\l£ 



2m 



1-5 



I/+ In 



ud 



VFq/'2 + l-k)B\ 



v^ In 



ud 



vpql2 - fioB\ 



(20) 



For q = and in the absence of the Zeeman field B the renormalized interaction diverges logarithmically. For finite q 
and B, the odd power terms in the expansion of £^/9+k have opposite signs, and hence, they cut off the logarithmic 
divergence. 

Because the density of states is larger in the — band, the interaction vertex is most divergent for: 



(clean) _ 2^oi3 
^ Vp 



(21) 



To linear order in the magnetic field and the pair momentum, this choice of q cancels the effect of the magnetic field 
on pairs in the — band. With this choice, the logarithmic divergence in the — branch is cut off only by terms of order 
O (B^/A^^). The critical field Be is the field beyond which the interaction vertex is finite: 



In Ao « iy+ Hf^oBc) + v^ In ^^^^ , 



(22) 



where Aq = ujp, exp [j^j^^^g 
From Eq. 22 we identify: 



u}p) exp \ ^ J is the gap in the absence of the exchange field, and A^o = 2akp. 



HqBc ~ Aq 



l + o'/fF 

A^\ 2+-/-. 
Ao 



(23) 



Anisotropic Bands 



We now re-examine the possibility of finite-momentum pairing in systems with anisotropic band dispersions. This 
is relevant, for example, in the oxide interface system where we expect superconductivity to arise predominantly in 



bands descending from d^z and dyz orbitals[? ], which have anisotropic hopping amphtudes in the xy plane. We 
model the band-anisotropy by introducing different effective masses for the x and y directions, replacing the kinetic 
energy term by: 

^^:^ + :^-^{l - ecos2^fc) (24) 

2m Zirix itJiy Zm 

where m = ^ ^^ , e — ^, ^ ,^ ^ e ("l? 1) is the eccentricity of the kinetic term, and (f>k is the angle between the 
momentum and the x-axis. 

The effect of anisotropic bands is easiest to analyze in the limits Ago 3> /i, and A^o <^ fi. For the case of Ago 3> /i, 
the problem is essentially identical to the isotropic case described above, once we replace m — > my in Eq. 21123 In 



this limit, the optimal pairing wave- vector becomes Oopt = ^^2_^ and the critical field is related to that of the isotropic 

band case by: hqB^™^- « {^) ^J■oB^c°' ■ In this limit, the critical field is either enhanced (suppressed) relative to 

the isotropic case if the band has larger (smaller) band-mass in the direction perpendicular to the exchange field, B. 

In the opposite limit, Ago <S^ /i, one can find a closed form expression for the critical field Be for weak anisotropy 

(e <C 1): i?^'"*^ « (1 -I- ^ j i3Jf°. In this limit, the critical field is enhanced by the anisotropy independent of the 

direction of the anisotropy. This can be understood simply by noting that the anisotropic dispersion leads to an 
elongated Fermi-surface, leading to an enhancement of nesting in the Cooper channel for pairs forming along the 
elongated direction. 

Generically, we find that the anisotropy can shift the quantitative values of the pairing vector gopt and the critical 
field Be, but do not alter the parametric enhancement of the critical field due to spin-orbit. We expect that band- 
anisotropy will be similarly unimportant in the presence of disorder. 

PAIRING INTERACTION VERTEX WITH DISORDER 

Having described pairing in the absence of disorder, we now examine how the analysis is modified in the presence of 
disorder. Scattering by impurities broadens the bands, giving the quasiparticles finite (elastic) life-times r. The finite 
life-time contains contributions from inter- as well as intra-band scattering. From the general theory of superconduc- 
tivity we know that disorder does not affect the superconducting state as long as the Cooper pairs remain correlated 
after each scattering event. In other words, if the probability of an electron in state k to be scattered into state k' 
is equal to the probability of its time reversed partner to be scattered into the partner of k'. Quantitatively, these 



correlation are described by the Cooperon (see Eq. 16 and subsequent discussion). Fully correlated pair scattering 
manifests itself in a massless Cooperon, i.e., the Cooperon, having a diffusive pole, diverges at zero frequency and 
momentum. As long as A^qT < 1, the pair scattering is dominated by processes where both electrons either remain in 
the same Rashba band or transfer to the other band, and, in principle, the Cooperon could have a diffusive pole. In 
this scenario, the superconducting state would not be sensitive to disorder. However, as we show below, the magnetic 
field reduces the pair correlation after certain scattering events, and the Cooperon is massive. This effect is most 
relevant for weak disorder. 

An additional effect of disorder is related to spin precession. As a consequence of the spin orbit coupling, in the 
absence of magnetic field, the spin of an electron with momentum k precesses with frequency Ago around the vector 
k. Since time reversal symmetry is preserved, electrons with ±k can pair to create a superconducting state. When 
a magnetic field is applied and time reversal symmetry is broken, the magnetic field changes the magnitude and 
direction of the spin precession (we consider here only the limit ij^qB <C Ago where the spin precession due to the 
spin-orbit coupling is weakly modified by the magnetic field). Correspondingly, the electrons with momenta ±k are 
no longer related by time reversal (the ±k pairs dephase), and the magnetic field serves as pair breaking mechanism. 
Note that changing the pair momentum ±k — >■ ±k -|- q/2 not only effect the spin precession, it also changes the kinetic 
energy of the pair. Superconductivity is obtained when these two dephasing mechanisms are minimized together. 
As a consequence of scattering by impurities the spin precession relaxes. Each time an electron is scattered into a 
different momentum state, it experiences a different torque which leads to relaxation of the precession. For disorder 
weaker than the spin orbit coupling, the large (time reversal preserving) component of the spin precession is almost 
not affected by the impurities. Nevertheless, the additional (breaking time reversal) precession of the spin due to the 
magnetic field can be suppressed by disorder. The latter, being relevant in the regime ^qB < t~^, actually helps the 
superconducting state to survive at high magnetic fields, as we show below. 

Taking all elements into account, we identify three distinct regimes: 1) weak disorder Aqt > 1, 2) moderate disorder 
Aq < T~^ < Ago, and 3) strong disorder A^qT < 1, and analyze each in turn. The first two regimes were studied in 
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[5] in the limit of weak spin-orbit coupling. Here we extend the analysis to all regimes to treat arbitrary strengths of 
spin-orbit and disorder. 

Weak Disorder (tAq > 1) 

For weak disorder, t~^ ^ t^ « IJ-oBc, we can assume that just like in the clean case, \vFq/'2 — fJ-oB] < t~^, 
while \vFq/2 -I- fJ,oB\ > t~^ (as we will see this assumption is consistent with our solution). Then, the pair coherence 
is strongly suppressed after scattering from the — to the + band and the Cooperon is massive (with mass of the 



order r ^). Consequently the Cooperon correction can be neglected, and Eq. 19 is only modified by the quasiparticle 



life-time entering the single particle Green's functions (see Eq. 17). Integrating over the frequency and momentum 



in Eq. 19 we find that F diverges when: 



-! = :._ In , "^^ = + v+ In , '^^ ^. (25) 

9 1/2t + v/1/4t2 + {vpq - 2fioB)^) 1/2t + ^1/4t^ + [vpq + 2^0^)2) 

The optimal pairing momentum is identical to the one obtained in the clean limit qopt — 2fiBy/vF- The critical field, 

on the other hand, is strongly suppressed, (J-qBc ~ Ao(Aot) i-°/"f . This expression for the critical field is valid as 
long as {fj,oBc)^/Al^ < 1/r < Be- The critical exchange field decreases as the disorder increases, and as the system 
approaches the limit 1/t « Aq the critical field gets closer to the Pauli limit, /io^c ~ Aq. 

In this weak disorder limit, the fact that the quasiparticle life-time enters the gap equation through the single 
particle Green's function for the — band, while the Cooperon is massive tells us that some scattering events suppress 
pair coherence. More precisely, as long as the two electrons forming a Cooper pair scatter between states in the — 
band, they remain correlated. However, the probability of an electron to be scattered from one band to the other is 
not negligible. Scattering of a Cooper pair from the — to the -I- band (where the magnetic field serves as a strong 
pair breaking mechanism) , results in decoherence of the pair. [One can check that if scattering is limited to a single 
Rashba band, the FFLO state is immune to weak disorder.] 

Moderate Disorder (Ao < r^^ < Aso) — "Diffusive Regime" 

In the previous section we saw that weak disorder strongly suppresses the critical field as r^^ — > Aq. One might 
expect that stronger disorder could only further weaken the FFLO state. Surprisingly the finite momentum pairing 
re-emerges for moderate disorder Aq < t^^ < Ago, and the critical field increases beyond Aq. This regime has been 
studied in Ref. [2[ for the weak spin orbit limit, a <ti vp. The key element is that tb > t. Then, the component 
of the spin precession caused by the magnetic field is diffusive, and hence, strongly suppressed. Similar to the effect 
of scattering on momentum, the frequent changes in torque after each scattering event leads to diffusion of the 
vector of spin precession. As a consequence of the diffusion of both spin precession and momentum, the mechanism 
that previously reduced significantly the pair coherence gets weaker with increasing disorder. Correspondingly, the 
magnetic field enters the Cooperon as r^ r <?C t~^, i.e. in the same way as the pair momentum q. Therefore, the 
effect of B is smaller than in the weak disorder limit, and the system can tolerate higher magnetic fields. This is the 
main difference between the FFLO state in systems with and without spin-orbit coupling. In the absence of spin-orbit 
coupling, the magnetic field appears in the Cooperon as 1/tb- 

Incorporating the Cooperon in the vertex equation, the instability point is determined by: 



9 ^ l-V^G+G+{{))/2-V^G-G-{Q)/2 + V^[G+G+{l)-G-G-{l)Y/2' 



In the above equation and below we use the notation G\G\i{m) for the m-th angular-harmonic of the two Green's 
function averaged over the momentum: 



^^e-™«GA(q/2 + ki, £„)Gv(q/2 - ki, -e 



GAGA'(m)- / ^^e-™^GA(q/2 + ki,£„)GA-(q/2-ki,-£„). (27) 



Performing the frequency integration in Eq. |26[ we get: 



In Ao = In 



^{vpq + 2^ioBfT + ^{vpq - 2^ioBfT + -{aq - 2/ioB)V 



(28) 
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The instability still occurs at finite momentum but the pair momentum in the diffusive regime is different than the 
one obtained in the clean case: 



(diff) _ 4:a^loB 
9opt - ^2 I „2 



(29) 



In this regime, the optimal q is not simply given by optimizing one Rashba branch at the expense of the other. Rather, 
it is chosen to compromise between both branches in order to minimize the frequency in which the Coopcron diverges. 
If both branches had the same densities of states, pairing would occur at zero momentum. In the limit a ^ w^?, we 
recover the result of Ref . [5] . 



The critical field corresponding to g^pt is Be 
ering its value in the clean limit as t^^ — > Aso- 



Ao(Aot) ^", which grows with increasing disorder almost recov- 



Strong Disorder (r ^ ^ Aso) 

The revival of the FFLO state in the limit of moderate disorder can be attributed to the diffusive behavior of 
both the electron momentum and spin. For even stronger impurity scattering 1/t 3> Aso, the precession picture 
described previously breaks down. More concretely, while the momentum relaxation remains to be r~^, the spin 
relaxation becomes smaller due to motional narrowing f~^^ « ^so''" (the Dyakonov-Perel mechanism [3]). In the 
extreme limit A^qT — >■ 0, the direction of the spins decouples from the momentum and the spin-orbit coupling is 
completely ineffective. We then expect the critical field to approach the Pauli limit. 

These expectations are confirmed in calculating the vertex of the attractive interaction. Indeed, in the strong 
disorder limit, the Cooperon loses its diffusive behavior, and the gap equation is, 



r-i = 1 - .g^[G+G+(0) + G_G_(0)] 



2y/{A-B)^-C^ 



A-B 



2v/(A-B)2-G2 
A^B~ V{A-B) 



G2 



A + B-^/{A-B)^-C^ 



(30) 



where 



A: 

B 

C 



V 



2 r 



2 

4 



G+G+(0) + G_G_(0) 



aq 



2A, 



2noB y 



G+G+(0) + G_G_(0) + G_G+(0) + G+G_(0) 



G+G+(1)-G_G_(1) 



aq + 2/io-B 

2AZ, 



(31) 



Integrating Eq. 30 over the frequency, the condition for the onset of pairing becomes: 



lnAo = 



V^-' 



2 - {aq - 2^oS)2/4 

T 



A2 



so 



VM 



{aq - 2/ioB)2/4 



In 



In 



A 



(aq + 2^JLaBf 



A'^ T 

'—^so ' 



v/Air2-(ag-2Moi?)V4 



ALr + ^"^ + ^;^^ + VALr2 - (aq ^ 2,oB)yA 



(32) 



In the above expression we have dropped terms of the order (vpqT)'^, (Bt)^ <^ 1. For q — our expression coincides 
with the one derived in Ref. [1]. We can see that all of the natural choices for pairing momentum: q = and 
q — ±2fj,QB/a, give rise to critical fields of the same order: Be ~ Ao\/A2^t/Ao « Ao(Aofso)~^/2 which eventually 
decays to the Pauli limit as disorder becomes infinitely strong. To determine which of the three possible values of 
q is the optimal pair momentum, one has to examine more carefully the numerical coefficients in the corresponding 
expressions for Be- Nevertheless, as 1/t, becomes large enough, we expect to have a first order phase transition from 
the FFLO state into a superconducting state with q = 0. 
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DISCUSSION 

To summarize, wc have presented the derivation of the critical field as a function of disorder for four different 
regimes. We have shown that scattering by impurities strongly affects the superconducting state with finite pairing 
momentum. While weak disorder significantly reduces the critical exchange field, the diffusive behavior obtained 
when 1/r > Aq helps to increase the field almost up to its value in the clean limit. However, crossing to the regime 
of strong disorder, 1/r > A^o, the direction of the spin is no longer correlated with the momentum and there is no 
difference between a system with and without spin-orbit coupling. Therefore, the FFLO state disappears at strong 
disorder, and the critical field approaches the Pauli limit. 
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